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Abstract 

Aims. We present the analysis of the gravitational microlensing event OGLE-2011-BLG-0251. This anomalous event was observed 
by several survey and follow-up collaborations conducting microlensing observations towards the Galactic Bulge. 
Methods. Based on detailed modelling of the observed light curve, we find that the lens is composed of two masses with a mass 
ratio q = 1.9 x 10"'^. Thanks to our detection of higher-order effects on the light curve due to the Earth's orbital motion and the 
finite size of source, we are able to measure the mass and distance to the lens unambiguously. 

Results. We find that the lens is made up of a planet of mass 0.53 ± 0.21 Mj orbiting an M dwarf host star with a mass of 
0.26 ± O.11M0. The planetary system is located at a distance of 2.57 ± 0.61 kpc towards the Galactic Centre. The projected 
separation of the planet from its host star is d = 1.408 ± 0.019, in units of the Einstein radius, which corresponds to 2.72 ± 0.75 
AU in physical units. We also identified a competitive model with similar planet and host star masses, but with a smaller orbital 
radius of 1.50 ± 0.50 AU. The planet is therefore located beyond the snow line of its host star, which we estimate to be around 
~ 1 - 1.5 AU. 
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1. Introduction 

Gravitational microlensing is one of the methods that al- 
lows us to probe the populations of extrasolar planets 
in the Milky Way, and has now led to the discoveries 
of 16 planet^, several of which cou ld not have been de- 
tected with other techniques (e.g. Beaulieu et alJ l2006l 
iGaudi et al] l2008l iMuraki et all |201lD. In particular, mi- 
crolensing events can reveal cool, low-mass planets that 
are difficult to detect with other methods. Although this 
method presents several observational and technical chal- 
lenges, it has recen t ly led to several significant scientific 
results. ISumi et aD (l201l[ ) analysed short time-scale mi- 
crolensing events and concluded that these events were pro- 
duced by a population of Jupitcr-mass free-floating planets, 
and were able t o estimate the numb er of such objects in 
the Milky Way. ICassan et al.l ()2012D used 6 years of ob- 
servational dat a from the PLANE T co llaboration to buil d 
on the work of iGould et all (|2010D and lSumi etHI (|201lD . 
and derived a cool planet mass function, suggesting that, 
on average, the number of planets per star is expected to 
be more than 1. 

Modelling gravitational microlensing events has been 
and remains a significant challenge, due to a complex 
parameter space and computationally demanding cal- 
culat ions. Recent developments in mod e lling methods 
(e.g. ICassanI [2008 1: 'Kains ct al.''200l [20T1 iBennettI [20Tot 
iRvu et al.ll2010l : [Bozza ct al 2012), however, have allowed 
microlensing observing campaigns to optimise their strate- 
gies and scientific output, thanks to real-time modelling 
providing prompt feedback to observers as to the possible 
nature of ongoing events. 

In this paper we present an analysis of microlensing 
event OGLE-2011-BLG-0251, an anomalous event discov- 
ered during the 2011 season by the OGLE collaboration 
and observed intensively by follow-up teams. In Sec. ^ we 
briefly summarise the basics of relevant microlensing for- 
malism, while we discuss our data and reduction in Sec. [31 
Our modelling approach and results are outlined in Sec. S) 
we translate this into physical parameters of the lens sys- 
tem in Sec. [S] and discuss the properties of the planetary 
system we infer. 



2. Microlensing formalism 

Microlensing can be observed when a source becomes suf- 
ficiently aligned with a lens along the line of sight that the 
deflection of the source light by the lens is significant. A 
characteristic separation at which this occurs is the Einstein 
ring radius. When a single point source approaches a sin- 
gle point lens of mass M with a projected source-lens sep- 
aration u, the source brightness is magnified following a 
symmetric "point source-point lens" (PSPL) pattern which 
can be parameterised with an impact parameter uq and a 
timesc ale ^f,i both exp ressed in units of the angular Einstein 
radius ()Einsteinlll936l ). 



4GM ZDs - 



(1) 



where G is the gravitational constant, c is the speed of light, 
and Dg and Dl are the distances to the source and the 
lens, respectively, from the observer. The timescale is then 
= ^e/Mj where fi is the lens-source relative proper mo- 
tion. Therefore the observable is a degenerate function 
of M, Z)l and the source's transverse velocity v± , assuming 
that Ds is known. However, measuring certain second-order 
effects in microlensing light curves such as the parallax due 
to the Earth's orbit allows us to break this degeneracy and 
therefore measure the properties of the Icnsing system di- 
rectly. 

When the lens is made up of two components, the mag- 
nification pattern can follow many different morphologies, 
because of singularities in the lens equation. These lead 
to source positions, along closed caustic curves, where the 
lensing magnification is formally infinite for point sources, 
although the finite size of sources means that, in prac- 
tice, the magnification gradient is large rather than infi- 
nite. A point-source binary-lens (PSBL) light curve is of- 
ten described by 6 parameters: the time at which the source 
passes closest to the center of mass of the binary lens, to, 
the Einstein radius crossing time, tE, the minimum impact 
parameter uq, which are also used to describe PSPL light 
curves, as well as the source's trajectory angle a with re- 
spect to the lens components, the separation between the 
two mass components, d, and their mass ratio q. Finite 
source size effects can be parameterised in a number of 
ways, usually by defining the angular size of the source 
in units of 6'e : 



(2) 
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where 9^ is the angular size of the source in standard units. 



3. Observational data 

The microlensing event OGLE-2011-BLG-0251 was 
discovered by the OGLE (Optical Gravitational Lens 
Experiment) collaboration's Early Warning System 
()Udalskil |2003|) as part of the release of the first 431 
microlensing alerts following the OGLE-IV upgrade. 
The source of the event has equatorial coordinates 
a = 17''38™14.18^ and S = -27°08'10.1" (J2000.0), or 
Galactic coordinates of 6) = (0. 670°, 2.334°). Anomalous 
behaviour was first detected and alerted on August 9, 2011 
(HJD~'2455782.5) thanks to real-time modelling efforts by 
various follow-up teams that were observing the event, but 
by that time a significant part of the anomaly had already 
passed, with sub-optimal coverage due to unfavourable 
weather conditions. The anomaly appears as a two-day 
feature spanning HJD = 2455779.5 to 2455781.5, just 
before the time of closest approach Iq. Despite difficult 
weather and moonlight conditions, the anomaly was 
securely covered by data from five follow-up telescopes in 
Brazil (//FUN Pico dos Bias), Chile (MiNDSTEp Danish 
1.54m) New Zealand (/iFUN Vintage Lane, and MOA Mt. 
John B&C), and the Canary Islands (RoboNct Liverpool 
Telescope). 

The descending part of the light curve also suffered from 
the bright Moon, with the source ^ 5 degrees from the 
Moon at ~ 85% of full illumination, leading to high back- 
ground counts in images and more scatter in the reduced 
data. We opted not to include data from Mt. Canopus Im 
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telescope in the modelling because of technical issues at the 
telescope affecting the reliability of the images, and also ex- 
cluded the /-band data from CTIO because they also suf- 
fer from large scatter, probably due to the proximity of the 
bright fuU Moon to the source. 

The data set amounts to 3738 images from 13 sites, from 
the OGLE survey team, the MiNDSTEp consortium, the 
RoboNet team, as weh as the //FUN, PLANET and MOA 
collaborations in the /, V and R bands, as well as some 
unfiltered data; data sets are summarised in Table [T] and 
the light curve is shown in Fig. [TJ We redu ced all data us- 
ing the difference im aging pipeline DanDIA ('Br amichll2008t 
iBramich et al.|[2"013[ ). except for the OGLE data, which was 
reduced by the OGLE team with their optimised offline 
pipeline. 

For each data set, we applied an error bar rescaling fac- 
tors a and b to normalise error bars with respect to our 
best-fit model (see SecS]), using the simple scaling relation 
(T- = a^/af+^ where cr- is the rescaled error bar of the i^^ 
data point and cr^ is the original error bar. The error bar 
rescaling factors for each data set is given in Table [TJ We 
did not exclude outliers from our data sets, unless we had 
reasons to believe that an outlier had its origin in a bad 
observation, or in issues with the data reduction pipeline. 



4. Modelling 

We modelled the light curve of the event using a Markov 
Chain Monte Carlo (MCMC) algorithm with adaptive step 
size. We first used the "standard" PSBL parameterisation 
in our modelling, whereby a binary-lens light curve can be 
described by 6 parameters: those given in Sec. [51 ignoring 
the second-order parameter described in that section. 
For all models and configurations we searched the parame- 
ter space for solutions with both a positive and a negative 
impact parameter uq. 

We started without including second-order effects of the 
source having a finite size or parallax due to the orbital mo- 
tion of Earth around the Sun, and then added these sep- 
arately in subsequent modelling runs by fitting the source 
size parameter p*, as defined in Sec. [U and the parallax 
parameters described below. Both effects led to a large de- 
crease in the statistic of the model (> 1000), which could 
not be explained only by the extra number of parameters. 

For the finite-source effect, we additionally considered 
the limb-darkening variation of the source star surface 
brightness by modelling the surface-brightness profile as 



Iip.x = /o,a[1 - q (1 - cosV')] , 



(3) 



where /q,^ is the brightness at the centre of the source, and 
tjj is the angle between a normal to the surface and the line 
of sight. We adopt the limb-darkening coefficients based 
on the source type determined from the dereddened colour 
and brightness (see Sec. 15. ip . The values of the adopted 
coefficients are cy = 0.073, cj = 0.624, cji = 0.542, based 
on the catalogue o f lClared((200?il) . 

Finally, in a third round of modelling, we included both 
the effects of parallax and finite source size ( "ESBL -I- par- 
allax"). Including these effects together led to a significant 
improvement of the fit, with Ax^ > 500 compared to the 
fits in which those effec ts were added separately. Computing 
the /-statistic (see e.g. lLuptonlll993l ) for this difference tells 
us that the probability of this difference occurring solely 



due to the number of degrees of freedom decreasing by 1 or 
2 is highly unlikely. Our best-fit ESBL + parallax model is 
shown in Fig. [TJ 

To model the effect of parallax, we use d the geocen - 
tric formalism (|Dominiklll998l: lAn et al.ll2002l: lGouldll2004[) . 
which has the advantage of allowing us to obtain a good 
estimate of toi^E and uq from a fit that does not include 
parallax. This formalism adds a further 2 parallax param- 
eters, tte^e and t:e,n, the components of the lens parallax 
vector -ke projected on the sky along the east and north 
equatorial coordinates, respectively. The amplitude of tt^; 
is then 



■KE = 



'E.N ■ 



(4) 



Measuring tt^; in addition to the source size allows us to 
break the degeneracy between the mass, distance and trans- 
verse velocity of the lens system that is seen in Eq. ([TJ) . This 
is because -ke also relates to the lens and source parallaxes 
TTi and TTg as 



7r_B = 



(5) 



Using this in Eq. ([TJ) allows us to solve for the mass of the 
lens. 

As an additional second-order effect, we also consider 
the orbital motion of the binary lens. Under the approxima- 
tion that the change rates of the binary separation and the 
rotation of the binary axis are uniform during the event, the 
orbital effect is taken into consideration with 2 additional 
parameters of d and d, which represent the rate of change of 
the binary separation and the source trajectory angle with 
respect to the binary axis, respectively. It is found that the 
improvement of fits by the orbital effect is negligible and 
thus our best-fit model is based on a static binary lens. 

Below we outline our modelling efforts that resulted in 
fits that were not competitive with our best-fit ESBL -(- 
parallax models, and which we therefore excluded in our 
light curve interpretation. 

4.1. Excluded models 
4.1.1. Xallarap 

We attempted to model the effects of so-called xallarap, or- 
bital motion of the source if it has companion (iGriest fc Hul 
I1992D . ModcUing this requires five additional parameters: 
the components of the xallarap vector, ^e,n and ^e,e, the 
orbital period P, inclination i and the phase angle ip of 
the source orbital motion. By definition, the magnitude of 
the xallarap vector is the semi-major axis of the source's 
orbital motion with respect to the centre of mass, as, nor- 
malised by the projected Einstein radius onto the source 
plane, te = DqOe, i.e. 



(6) 



The value of as is then related to the semi-major axis of 
the binary by 



as = 



aM2 



(7) 



Ml + M2 ' 

where Mi and M2 are the masses of the source components 
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Figure 1. Light curve of OGLE-2011-BLG-0251. Data points arc plotted with 1-a error bars, and the upper panel shows 
a zoom around the perturbation region near the peak 



In Fig. [21 we show of the fit plotted as a function 
of the orbital period of the source star. We compare this 
to the statistic of the best parallax fit. We find that 
xallarap models provide fits competitive with the parallax 
planetary models for orbital periods P > 1 year. However, 
the solutions in this range cannot meet the constraint pro- 
vided by the source brightness. Combining Equations (jS]) 
and (E) with Kepl er's third law, 
(|Dong et all 120090 



upper limit for its mass of 3 Mq. With this constraint, we 
find that xallarap models are not competitive with parallax 
planetary models, and we therefore exclude the xallarap 
interpretation of the light curve. 



a7(Afi +M2) yields 412. Binary source 



AU 



Mi 



(8) 



Rearranging this equation for M2, and using the fact that 
M2/{Mi + M2) < 1, we can derive an upper limit for the 
mass of M2, 



2,min 



p2 



(9) 



In the lower panel of Fig. [21 we show the minimum mass 
of the source companion as a function of orbital period. 
The blending constraint means that the source companion 
cannot be arbitrarily massive, and we use a conservative 



We also attempted to model this event as a binary source 
- point lens (BSPL) event. For this we introduced three 
additional parameters: the impact parameter of the sec- 
ondary source component, uo,2, and its time of closest ap- 
proach, io,2, as well the flux ratio between the source com- 
ponents. We note that parallax is also considered in our 
binary source modelling, for fair comparison to other mod- 
els. We find that the best binary-source model provides a 
poorer fit, with = 3809, which gives Ax^ ~ 180 com- 
pared to our best planetary model (including parallax, see 
model D in the following section). Residuals for this model, 
as well as all other models discussed in this section are 
shown in Fig. [31 
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Table 1. Data sets for OGLE-2011-BLG-0251, with the number of data points for each telescope/ filter combination. 
The rcscaling coefficients a and b are also given, with error bars rescaled as a' = asj + 5^, where a' is the rescaled 
error bar and a is the original error bar. 
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Figure 2. Constraints from the xallarap fit as a function 
of the orbital period P of the source star. The top panel 
shows of the xallarap fit as a function of P, with a red 
circle marking the location of the best parallax model. The 
bottom panel shows the minimum mass of the source com- 
panion as a function of P. The shaded area in both panels 
indicates where models are excluded based on conservative 
blending constraints on the source companion's mass. 



the degeneracy among these local solutions is rather severe, 
as can be seen from the residuals shown in Fig. [31 

For the identified local minima, we then further refine 
the lensing parameters by conducting additional modelling, 
considering higher-order efi'ects of the finite source size and 
the Earths orbital motion. It is found that the higher-order 
effects are clearly detected with Ax^ > 500. Best-fit pa- 
rameter for each of the local minima are given in Table [21 
while Fig. [51 shows the geometry of the source trajecto- 
ries with respect to the caustics for all four minima. We 
note that the pairs of solutions A and D, and B and C, 
are degenerate under the well-known d o degeneracy 
([Griest fc Safizadeh|[T99llDominiklll999D : this is caused by 
the symmetry of the lens mapping between binaries with d 
and d~^. Comparing the pairs of solutions, we find that the 
A-D pair is favoured, with Ax^ > 40 compared to the B-C 
pair. On the other hand, the degeneracy between the A and 
D solutions is very severe, with only Ax^ ~ 7. In Fig. [HI we 
also show parameter-parameter correlations plots for model 
D, showing also the uncertainties in the measured lensing 
parameters. 



5. Lens Properties 

In this section we determine the properties of the lens sys- 
tem, using our best-fit model parameters, i.e. our wide- 
configuration ESBL -I- parallax model. We also calculated 
the lens properties for the competitive close-configuration 
model, with both sets of parameter values listed in Table [31 



4.2. Best-fit models 

We searched the parameter space using an MCMC algo- 
rithm as well as a grid of (d,q, a) to locate good st arting 
points for the algorithm (see e.g. iKains et al.ll2009h . over 
the range -4 < logg < and -1.0 < logd < 2. This 
encompasses both planetary and binary companions that 
might cause the central perturbation. In Fig. [31 we present 
the distribution in the d, q plane. We find four local so- 
lutions, all of which have a mass ratio corresponding to a 
planetary companion. We designate them as A, B, C and D; 



5.1. Source star and Einstein radius 

We determined the Einstein radius by first calculating the 
angular size of the source. This can be done by using th e 
magnitude and colour of the source (e.g. lYoo et all[200l . 
and empirical relations between these quantities and the 
angular source size. We start by using the location of the 
red giant clump (hereafter RC) on our colour-magnitude 
diagram (Fig. [7]) to estimate the reddening and extinction 
along the line of sight. We use an /-band abs olute magni- 
tude for the RC of M/,rc,o -0.12 ± 0.09 (iNataf et all 
I2012D . as weU as a colour {V - I)nc.o = 1-06 ± 0.12 
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d 


0.775 ± 0.010 


0.997 ± 0.009 


1.066 ± 0.001 


1.408 ± 0.019 




(1.68 ± 0.11) 


(0.93 ± 0.03) 


(1.11 ± 0.06) 


(1.92 ± 0.12) 




-0.33 ± 0.04 


-0.37 ± 0.04 


-0.40 ± 0.05 


-0.34 ± 0.05 




0.09 ± 0.01 


0.08 ± 0.01 


0.08 ± 0.01 


0.09 ± 0.01 




0.34 ± 0.04 


0.38 ± 0.04 


0.41 ±0.05 


0.35 ± 0.05 


"ff = Fb/Fs 


0.387 ± 0.035 


0.394 ± 0.001 


0.394 ± 0.042 


0.376 ± 0.017 


"is 


15.99 ± 0.03 


15.98 ± 0.01 


15.98 ± 0.02 


15.97 ± 0.01 


a ra 


16.97 ± 0.07 


16.99 ± 0.01 


16.99 ± 0.06 


17.04 ± 0.03 



Table 2. Best- fit model parameters and l-cr error bars for the four identified best binary- lens models including the effects 
of the orbital motion of the Earth (parallax). MHJD=HJD-2450000. ''for the OGLE data set 



(|Bensbv et al.ll20lH) . We compare these values to those on 
our colour- magnitude diagram (CMD), which we generated 
using OGLE /— and V— band photometry. From Fig. [71 
the location of the RG on our CMD is 

(/, V - /)rc = (17.19 ± 0.05, 3.45 ± 0.05) (10) 

so, using a distance modulus oi fj, = 14.52 z fc 0.09, i.e. a 
dista nce to the Galactic bulge of 8.0 ± 0.3 kpc (jYelda et al.l 
[201I . we find Ai = 2.79±0.10 and £;(!/-/) = 2.39±0.13. 

Using these values, the best-fit value for the magnitude 
of the source Is = 15.97±0.01, a source colour (V — I)s.o = 
1.15, and the empirical relations of iKervella &: Fouaug 
(|2008D . we find an angular source radius 9^, = 10.41 ± 



1.18 ^as, or a source star radius of i?* = 10.53 ± 1.19 i?Q. 
This, together with the best-fit value of the source size pa- 
rameter p*, allows us to calculate the size of the Einstein 
radius, 9^ — 9^,1 p^,. Using the relevant parameter val- 
ues, we find 6Ie = 0.749 ± 0.283 mas. This in turn al- 
lows us to calculate the source-lens relative proper motion, 
= 9^/1^ = 4.28 ± 1.62 mas/yr. 

5.2. Masses of the Lens Components 

Gombining Eq. ([T]) and Eq. ([5]) allows us to derive an ex- 
pression for the mass as a function of the parallax vector 
magnitude tte (defined by Eq. H]) : 
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1 


-1 
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-3.5 - 
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-0.2 -0.1 0.1 

log d 
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Figure 4. map in the d, q plane, showing the location of the four local minima identified by our modelling runs. Out 
of these, local minima A and D are competitive, with local minima B and C having Ax^ ~ 50 and 70 respectively, for 
the same number of parameters. Minima A and D correspond to the close and wide ESBL + parallax models discussed in 
the text. Different colours correspond to Ax^ < 25 (red), 100 (yellow), 225 (green), and 400 (blue); we note that the 
map is based on the original data, before error-bar normalisation, and therefore the Ax^ levels arc slightly different from 
those given in Table [2l The top panel shows the breadth of our parameter space exploration, encompassing planetary 
and non-planteray mass-ratio regimes, while the bottom panel shows a zoom on the region where our local minima are 
located. 



Mr 



(11) 



Using values found in the previous section, and our best- 
fit parallax parameter value t^e = 0.35 ± 0.05 yields a total 
lens mass Ml = 0.26 ± 0.10 M0. Using the best-fit mass 
ratio parameter value oi q = (1.92 ± 0.12) x 10^'^ yields 
component masses of 0.26 ± 0.11 M© and 0.53 ± 0.21 Mj, 
where Mj is the mass of Jupiter. 

5.3. Distance to the Lens 

We can also rearrange Eq. ([1]) to derive an expression for 
the distance to the lens Dl, 



A. = 



1 



AGAI 



(12) 



Using our parameter values as well as the lens mass derived 
thanks to our parallax measurement, we find a distance to 
the lens of -Dl = 2.57±0.61 kpc. This distance allows us to 
carry out a sanity check of the lens mass we derived in the 
previous section. By assuming that the contribution from 
the blended light comes from the lens, we can derive an 
upper limit to the /—band lens magnitude Mj using our 
best-fit blending parameter: 



Mi.L = nii,b - 5 logio -Dl - 10 - : 



(13) 
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Figure 5. Source trajectory geometry with respect to the caustics for all four local minima identified in Fig.lU the source 
size is marked as a red circle. 



where to/, 6 is the apparent /—band magnitude of the blend, 
Aj^i^ is the extinction between the observer and the lens, 
and Dl is in kpc. In practice, ^ since the lens 

is in front of the source, so we use the extreme scenario 
where Ajj^ ~ Aj to derive an upper brightness limit (lower 
limit on the magnitude) for the lens. We find this to be 
Af/,L = 2.19 ±0.53 mag, which corresponds to a maximum 
mass of the lens of M^^max = 1-65 ± O.23M0, assuming a 
main sequence star mass-luminosity relation, and assuming 
that the secondary lens component (i.e. the planet) does not 
contribute to the blended light. This is much larger than the 
value we derived in Scc. l5.2l for the mass of the primary lens 
component, which suggests that some blending comes from 
stars near the source rather than from the lens, although it 
is difficult to quantify this without an estimate of 

Finally, we can also use the distance to the lens and the 
size of the Einstein ring radius to calculate the projected 
separation r± between the lens components in AU. Using 
our best-fit projected angular separation d = 1.408 ±0.019, 
we find a projected (i.e. minimum) orbital radius r±^ = 
2.72 ±0.75 AU. 

We can compare this to an estimate of the location 
of the "snow line", which is the location at which water 
sublimated in the midplane of the protoplanetary disk, i.e. 
the distance at which the midplane had a temperature of 
/mid = 170 K (although oth er studies have noted tha t this 
temperature varies; see e.g. iPodolak fc Zucked 1200411 . The 
core accretion model of planet formation predicts that gi- 
ant planets form much more easily beyond the snow line, 
thanks to easier condensation of icy material and there- 
fore easier formation of large solid cor es in the early stages 
of the circumstellar disk's evolution. [Kennedy fc KenvonI 
(|2008|) modelled the evolution of the snow line's location, 
taking into account heating of the disk via accretion, as 
well as the influence of pre-main sequence stellar evolution. 
Using a rough extrapolation of their results, we estimate 
that the snow line (at t = IMyr) for the planetary host 
star in OGLE-2011-BLG-0251 is located at around - 1-1.5 
AU. We therefore conclude that OGLE-2011-BLG-0251Lb 
is a giant planet located beyond the snow line, with both of 
our competitive best-fit models yielding projected orbital 
radii larger than 1.5 AU. 





close 


wide 


Of [/lasj 


10.29 ±1.17 


10.41 ±1.18 


[mas] 


0.71 ±0.26 


0.75 ±0.28 


/irci [mas yr"^] 


4.09 ±1.50 


4.28 ±1.62 


Ml [Mq\ 


0.26 ±0.10 


0.26 ±0.11 


M2 [Mj] 


0.45 ±0.18 


0.53 ±0.21 




1.71 ±0.23 


1.65 ±0.23 


Dl [kpc] 


2.71 ±0.61 


2.57 ±0.61 


r± [AU] 


1.50 ±0.50 


2.72 ±0.75 



Table 3. Lens properties derived as detailed in Sec. [5l for 
both competitive parallax models. 



Wc list all the lens properties in Table |31 both for 
the best-fit model parameters that we have used above, 
and for the close-configuration model, for comparison. Lens 
properties derived using the close-configuration model are 
very similar to those we found using the wide-configuration 
model, the only major difference being in the orbital radius. 
For the close model, we find an orbital radius of 1.50 ±0.50 
AU; which is close to the location of the snow line. 

6. Conclusions 

Our coverage and analysis of OGLE-2011-BLG-0251 has 
allowed us to locate and constrain a best-fit binary-lens 
model corresponding to an M star being orbited by a gi- 
ant planet. This was possible through a broad exploration 
of the parameters both in real time, thanks to the recent 
developments in microlensing modelling algorithms, and 
after the source had returned to its baseline magnitude. 
Various second-order effects, as well as other possible, non- 
planetary, interpretations for the anomaly were considered 
during the modelling process. Based on the best-fit solu- 
tion, we were able to constrain the masses and separation 
of the lens components, as well as various other character- 
istics, thanks to a strong detection of parallax effects due 
to the Earth's orbit around the Sun, in conjunction with 
the detection of finite source size effects. We found a planet 
of mass 0.53 ± 0.21 Mj orbiting a lens of 0.26 ± O.llAfg 
at a projected radius r± = 2.72 ± 0.75 AU; the whole 
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Figure 6. Parameter-parameter correlations for our 9 fitted parameters. Colours indicate the limits of the 1, 2, 3, 4 and 
5-cr confidence limits for each pairwise distribution. A closer view of the correlation between parallax parameters is shown 
on the top right inset. 



system is located at a distance of 2.57 ± 0.61 kpc. Our 
competitive second-best model leads to similar properties, 
but a smaller projected orbital radius = 1.50 ± 0.50. 
The two best-fit models are competitive and therefore we 
cannot make a strong claim about which orbital radius is 
favoured. However, comparing both values of the projected 
orbital radius to the approximate location of the snow line 
for a typical star of the mass of the primary lens compo- 
nent, we conclude that OGLE-2011-BLG-0251Lb is a gi- 
ant planet located around or beyond the snow line. This 
is in line with predictions from the core accretion model 
of planet formation, from which we expect large planets to 
be more numerous beyond the snow line; this is also where 
microlensing detection sensitivity is at its highest, enabling 
us to probe a region of planetary parameter space that is 
difficult to reach for other methods. 
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